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The prediction of long-term geothermal reservoir perfor- 
mance and the environmental impact of exploiting this resource are 
two important problems associated with the utilization of geothermal 
energy for power production. Our research effort addresses these 
problems through numerical simulation. Computer codes based on 
the solution of partial-differential equations using finite -element 
techniques are being prepared to simulate multiphase energy 
transport, energy transport in fractured porous reservoirs, well 
bore phenomena, and subsidence. 


I. INTRODUCTION 

Geothermal reservoir simulation, in general, refers to reproducing and 
predicting the behavior of a geothermal reservoir by use of a model. In this 
discussion we consider only mathematical models which we define as: (1) a set 
of equations that describe the physical processes active in a geothermal reser- 
voir, and (2) the solution of these equations subject to boundary and initial 
conditions. Because of the complex nature of the equations and boundary con- 
ditions, the solution generally involves numerical techniques in conjunction 
with the digital computer. 

One simulates a geothermal reservoir to estimate the quantity of recov- 
erable energy and the rate at which mass and energy may be extracted. To 
achieve these objectives, known geological information obtained from both 
surface techniques and drilling is utilized to determine equation parameters, 
boundary, and initial conditions. During the earlier stages of field develop- 
ment, this information is limited, and the simulation model will undoubtedly be 
crude. Through the simulation of a variety of producing schemes, however, 
the energy available and its rate of extraction may be roughly estimated. If the' 
decision is made to develop the field further, the reservoir model may be used 
to help answer such engineering questions as optimal well location and spacing, 
or whether or not to reinject condensate. As more geological information 
becomes available through continued drilling, the simulation model may be 
updated to give a more accurate analysis of the geothermal reservoir. 
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A. Previous Investigation . * 


The literature on the simulation of both heat and fluid flow in porous 
media is extensive. Much of the literature is concerned with the extraction of 
petroleum by in situ combustion or hot -fluid injection (e.g.. Refs. 1 and 2), 
Another area of interest is that associated with free convection models (see, 
for example, Refs. 3 and 4) . If, however, only geothermal reservoir simu- 
lation is considered, the literature is relatively limited. 

Whiting and Ramey (Ref. 5) were the first to apply a reservoir model to a 
geothermal system. Their model allowed for two-phase (steam-water) flow, 
but did not take into consideration the spatial dependence of the solution (it was 
a lumped parameter model), Harlow and Pracht (Ref. 6) developed a single 
phase (hot-water) model coupled with a rock-fracture model. They attempted 
to demonstrate that geothermal energy could be extracted from hot dry rocks. 
Mercer (Ref. 7) applied a hot-water model to the Wairakei hydrothermal 
system. The model was limited spatially to the two horizontal dimensions and 
was unable to reproduce historical data past 1962 due to the formation of steam 
in the reservoir. Brigham and Morrow (Ref. 8), in an attempt to allow some 
spatial variation in their solution, developed three lumped-parameter models 
based on vapor liquid distribution. Toronyi (Ref. 9) developed a two-dimensional 
(areal or cross-section) reservoir model coupled with a wellbore model. His 
reservoir model allowed only two -phase flow, being restricted to the saturated- 
vapor pressure curve. 

B. A General Simulator 


The above models could be thought of as contributions toward a general 
model capable of completely describing the response of a reservoir to exploita- 
tion. The development of such a model constitutes the principal objective of 
our current research effort. 

The physical properties of the fluid in a geothermal reservoir vary spa- 
tially and temporally to the extent that vertical velocities are often of the same 
order of magnitude as horizontal velocities. Consequently, a general model 
must be a distributed parameter model with the capability of three-dimensional 
simulation. 

To answer important engineering questions, a general model must not 
only accommodate the transient flow of compressed hot water, steam-water 
mixtures, and super-heated steam, but also must allow for phase changes 
(see Fig. 1). Presently we are neglecting species transport and chemical 
reactions that may occur in the geothermal system. 

Subsidence of the land surface is often encountered when large quantities 
of fluid are removed from subsurface reservoirs. Consequently, a general 
model should have the capability of forecasting the spatial and temporal dis- 
tribution of subsidence resulting from geothermal development. Such a simu- 
lator must solve not only the fluid flow and energy transport equations but also 
the equations describing the elastic behavior of the reservoir skeleton. 
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The governing equations employed in reservoir simulation are generally 
based on the physics of flow through porous media. While this approach may 
hold in selected cases, there are many reservoirs that are productive only 
because large fractures are encountered during drilling. To simulate such a 
fractured reservoir, the classical continuum approach to porous media flow 
must be modified and extended to accommodate the influence of fracture flow. 

While we are currently developing models which will include each of the 
important capabilities described above, we cannot herein describe all aspects 
of our research program. Consequently, we have elected to focus our attention 
on multi-phase energy transport, a topic which should be of interest to anyone 
contemplating the development of a geothermal reservoir simulation capability. 


II. THEORETICAL DEVELOPMENT 

The general governing equations consist of mass, momentum, and energy 
balances for each phase present in a geothermal system. 

A. Mass Balance 

The mass balances for steam, s, and water, w, may be written as: 
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v = phase average velocity, [Lt”^J 
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P = average density, [ML' ] 
q = source term, [ML - t” ] 
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d v = rate of vaporization, [ML" t - ] 

<(> = porosity, dimensionless 
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B. Momentum Balance 


For velocity, we assume that Darcy's equation for multiphase flow may 
be used (see Ref. 10) : 
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where 


k = local intrinsic permeability tensor, [L ] 
k y = relative permeability, dimensionless 
p = dynamic viscosity, [ML *t 
p = pressure, [ML" t ] 

-2 

g = gravitational acceleration, [Lt ] 
Combining Eqs. (1) with (3) and (2) with (4) gives 
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These equations describe the flow of steam and water in a porous medium 
where a pressure difference exists between the two phases. This pressure 
difference is generally defined as capillary pressure, p c> where 
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Probably the most important effect that capillary pressure has in a geothermal 
reservoir is to lower the vapor-pressure curve. Ramey, et al. (Ref. 11) point 
out that the reason for the lowering of the vapor-pressure curve is that vapor- 
pressure data found in steam tables (Refs. 12 and 13) are based on flat steam- 
water interfaces, whereas capillary pressure causes the interface in porous 
media to be curved. The amount the vapor-pressure curve is lowered in a 
geothermal reservoir is not completely understood. The work of Calhoun, 
et al. (Ref. 14), on consolidated rock does show a lowering of the vapor-pressure 
curve with decreased fluid saturation. The efforts of Cady (Ref. 15) and 
Bilhartz (Ref. 16), however, indicate no significant vapor pressure lowering in 
their experiments using unconsolidated sands. An important difference in the 
two results is that the experiments of Calhoun, et al. , were made at a temper- 
ature of 36° C, and those conducted by Cady and Bilhartz ranged from approxi- 
mately 121 to 240° C. Further work on the importance of capillary pressure in 
geothermal reservoirs is required. In this paper, capillary pressure is 
assumed to be negligible. Under this assumption Eqs, (5) and (6) combine to 
reduce the number of equations and unknown parameters by one: 
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where p is the density of the total steam water mixture, defined as: 


p = p S + p S (9) 

r s s r w w ' ' 


Note that in combining Eqs. (5) and (6) the vaporization term, d , has also been 
eliminated. v 

C. Energy Balance 


Three energy balances, analogous to the mass balances written for steam 
and water, can be written for steam, water, and rock. If it is assumed that the 
movement of fluid (water and steam) through porous media is sufficiently slow, 
and the surface areas of all phases concerned are sufficiently large, then thermal 
equilibrium exists between the rock, steam, and water. This assumption 
reduces the three energy-balance equations to the following one: 
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The term described by Eq. (11) represents the compressible work term, 
and is often assumed negligible and omitted from the energy balance. 
Other definitions of terms in Eq. (10) include: 

h = enthalpy of saturated steam, [L t” ] 
s 

h^ = enthalpy of saturated water, [L t~ ] 

T = temperature, [T] 

= -3 -li 

K = thermal dispersion tensor for the medium' [MLt T J 
m 

- 3 

= average rock density, [ML - ] 
h^ = rock enthalpy, [L t~ ] 

h = total enthalpy of the mixture, as defined by 
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D, Assumptions 

By making the assumptions of thermal equilibrium and negligible capillary 
pressure, we are left with two equations. We must obtain additional relation- 
ships for the remaining dependent variables as functions of these two. 

The two dependent variables specified as unknowns must uniquely define 
the thermodynamic state of the system. Further, it is desirable that these 
variables be commonly measured in a field situation. For these reasons, we 



choose to put Eqs. (8) and (10) in terms of the dependent variables h, the total 
enthalpy of the mixture, and p, the fluid pressure. To do this we make the 
following assumptions: 

(1) Porosity, <j>, is a function of pressure alone. 

(2) Rock density, p^, may be treated as constant. 

(3) Total density, p, steam density, p s , and water density, p , are 

considered functions of pressure and enthalpy. * 

(4) Steam enthalpy, h s , and water enthalpy, h w> can be treated as 
functions of pressure. 

(5) Temperature, T, can be treated as a function of pressure in the 
two-phase region, and can be treated as a function of pressure and 
enthalpy in the single-phase region. 

(6) Rock enthalpy, h r , can be considered as a function of temperature. 

(7) Viscosities, p. s andp. w , are considered functions of -temperature . 

(8) Saturations, S w and S s , are functions of enthalpy and pressure. 

(9) Relative permeability can be treated as a function of saturation 
alone and relationships similar to those in Brooks and Corey 
(Ref. 17) may be used. 

(10) Thermal dispersion tensor can be generalized and treated as a 
property of the medium. 

Ramey, et al. (Ref. 18), point out that relative permeability can also be 
a function of temperature. Using unconsolidated sand and working with oil and 
water, Poston, et al. (Ref. 19), observed that for increased temperatures, the 
relative permeability curves shift to the left on the saturation axis. Regarding 
the thermal dispersion tensor it should be noted that Mercer (Ref. 7) separates 
the medium thermal dispersion tensor into three parts: conduction in the solid 
phase, diffusion in the liquid phase, and a velocity-dependent dispersion in the 
liquid phase. Further, Somerton, et al. (Ref. 20), point out that the thermal 
conductivity of the medium is a function of temperature, porosity, and water 
saturation. 
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E. Final Equations 

Making the above assumptions, Eqs . (8) and (10) become: 



(14) 
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Equations (13) and (14) are in terms of the dependent variables enthalpy 
and pressure. As can be seen from Fig. 1, a pressure and enthalpy value for 
a given point in space will determine whether compressed hot water, a steam- 
water mixture or super-heated steam exists at that point. Thus, Eqs. (13) and 
(14) may be used to describe most types of geothermal reservoirs. To simulate 
the behavior of such a reservoir, however, a solution for Eqs. (13) and (14) is 
required. 

Equations (13) and (14) are nonlinear, and their solution requires the use 
of numerical techniques. Such techniques are beyond the scope of this paper 
and will be mentioned only briefly. They include finite difference techniques, 
the finite element method, or some combination of the two. The nonlinear 
coefficients are generally handled by iteration, using, for example, the Newton- 
Raphson method (Refs. 2 and 9). Finally, some form of weighting of the spa- 
tially dependent variables is commonly used, upstream weighting of perme- 
ability, for example (Ref. 2). 

III. CONCLUSIONS 

We have presented the governing equations for energy transport and the 
flow of water and steam in porous media in order to provide a basis for the 
mathematical modeling of geothermal reservoirs. The primitive form of these 
equations is straightforward and available in various forms in the literature. 

Our objective here was to take the general mass and energy balances and 
reduce them to workable equations. It should be emphasized that the method 
described represents only one of many possible approaches. The develop- 
ment of a method for obtaining solutions to these equations is a challenging 
problem and beyond the scope of this paper. 

This particular approach to geothermal modeling has been guided by 
several objectives: (1) the model must provide efficient and accurate solutions 
consistent with the available geological and hydrological data, (2) it should be 
applicable to hot-water, steam-water, and super-heated steam geothermal 
systems, and (3) the method must be easily applied to natural geothermal 
reservoirs. Two assumptions, peculiar to this development, were made to 
achieve these objectives: (1) thermal equilibrium exists between all phases, 
in steam, water, and rock; and (2) capillary pressure between the steam and 
water phases is negligible. Invoking these assumptions the basic relationships 
reduce to two partial differential equations written in terms of convenient field 
variables, fluid pressure and enthalpy. These final two equations in conjunc- 
tion with well-known thermodynamic relationships provide the basis for a 
numerical model which we feel fulfills our initial objectives. 

The suitability of this model for simulation of natural geothermal reser- 
voirs is limited in two general ways: (1) by the validity of the basic assump- 
tions (negligible capillary pressure, porous media, etc. ), and (2) by the avail- 
ability of reliable geological and hydrological and thermodynamical data for the 
reservoir to be simulated. Within these constraints, the suggested approach 
should prove useful in geothermal reservoir development and management. 
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